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ABSTRACT 

Radiatively inefficient accretion flows onto black holes are unstable due to 
both an outwardly decreasing entropy ( "convection" ) and an outwardly decreas- 
ing rotation rate (the "magnetorotational instability"; MRI). Using a linear mag- 
netohydrodynamic stability analysis, we show that long-wavelength modes with 
\/H 3> are primarily destabilized by the entropy gradient and that such 

"connective" modes transport angular momentum inwards (A is the wavelength of 
the mode, H is the scale height of the disk, and (3 is the ratio of the gas pressure to 
the magnetic pressure). Moreover, the stability criteria for the convective modes 
are the standard H0iland criteria of hydrodynamics. By contrast, shorter wave- 
length modes with X/H ~ /3 _1//2 are primarily destabilized by magnetic tension 
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and differential rotation. These "MRI" modes transport angular momentum out- 
wards. The convection-dominated accretion flow (CDAF) model, which has been 
proposed for radiatively inefficient accretion onto a black hole, posits that inward 
angular momentum transport and outward energy transport by long-wavelength 
convective fluctuations are crucial for determining the structure of the accretion 
flow. Our analysis suggests that the CDAF model is applicable to a magnetohy- 
drodynamic accretion flow provided the magnetic field saturates at a sufficiently 
sub-equipartition value (f3 >> 1), so that long- wavelength convective fluctuations 
with X/H 3> can fit inside the accretion disk. Numerical magnetohydro- 

dynamic simulations are required to determine whether such a sub-equipartition 
field is in fact obtained. 

Subject Headings: accretion — accretion disks — black hole physics — convection 
— instabilities — MHD — turbulence 



1. Introduction 

Models of radiatively inefficient accretion flows provide a useful framework for interpret- 
ing observations of low-luminosity black hole X-ray binaries and active galactic nuclei (see, 
e.g., Narayan, Mahadevan & Quataert 1998; Quataert 2001; Narayan 2002 for reviews). In 
the past few years, there has been rapid theoretical progress in understanding the dynamics 
of such flows. Much of this advance has been driven by numerical simulations. In particular, 
a number of hydrodynamic simulations have been reported in which angular momentum 
transport is put in "by hand" using an a prescription (e.g., two-dimensional simulations by 
Igumenshchev, Chen, & Abramowicz 1996; Igumenshchev & Abramowicz 1999, 2000; Stone, 
Pringle, & Begelman 1999; and three-dimensional simulations by Igumenshchev, Abramowicz 
& Narayan 2000). Such simulations reveal a flow structure very different from advection- 
dominated accretion flow models (ADAFs) that were proposed to describe the structure of 
radiatively inefficient accretion flows (Ichimaru 1977; Narayan & Yi 1994; Abramowicz et 
al. 1995). In an ADAF, the gas accretes rapidly, and the gravitational potential energy of 
the accreting gas is stored as thermal energy and advected into the central black hole. By 
contrast, in the hydrodynamic simulations, the rate of mass accretion is much smaller and 
strong radial convection efficiently transports energy outwards. The simulations have been 
interpreted in terms of a "convection-dominated accretion flow" model (CDAF; Narayan, 
Igumenshchev & Abramowicz 2000; Quataert & Gruzinov 2000; Abramowicz et al. 2002). 
In a CDAF, convection simultaneously transports energy outwards and angular momentum 
inwards, strongly suppressing the accretion rate onto the central black hole. 
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The relevance of the hydrodynamic simulations, and the CDAF model derived from 
them, is unclear. This is primarily because of the ad hoc treatment of angular momentum 
transport. It is believed that angular momentum transport in accretion flows is primarily due 
to MHD turbulence initiated by the magnetorotational instability (MRI; Balbus & Hawley 
1991; BH91). Balbus & Hawley have argued that, because of the fundamental role played by 
magnetic fields, hydrodynamic models typically cannot describe the structure of the accretion 
flow (or differentially rotating systems more generally; e.g., Balbus & Hawley 1998; Balbus 
2000; 2001). They have applied this criticism in detail to the CDAF model (Hawley, Balbus, 
& Stone 2001; Balbus & Hawley 2002; hereafter BH02). In an independent argument, BH02 
also suggest that CDAF models violate the second law of thermodynamics. We do not 
consider this issue here, but intend to deal with it in a separate paper. 

In this paper we discuss the physics of the CDAF model within the framework of MHD. 
Apart from clarifying the theoretical underpinnings of the CDAF model, we believe that our 
analysis also provides a useful illustration of the conditions under which a hydrodynamic 
model is applicable to an intrinsically magnetohydrodynamic situation; this is important in 
other astrophysical contexts, e.g., the excitation of density waves at Lindblad resonances in 
a magnetized accretion flow, or diskoseismology models of QPOs. 

The paper is organized as follows. In §2 we discuss the linear stability of a differentially 
rotating and thermally stratified plasma in the presence of a weak vertical magnetic field. 
We identify and explain the difference between "convective" modes and "MRI" modes and 
show that the former transport angular momentum inwards while the latter usually transport 
angular momentum outwards. Then, in §2.1 we generalize this result to an arbitrary field 
orientation and arbitrary stratification. In §3 we discuss the implications of our results. 
In particular, we use the linear analysis to speculate about the conditions under which the 
saturated nonlinear turbulence in an accretion flow might give rise to a CDAF-like structure. 
In §4 we summarize our results. 



2. Convective Modes in a Magnetized Differentially Rotating Plasma 

We consider the linear stability of a differentially rotating and thermally stratified 
plasma in the presence of a weak magnetic field. We specialize to the case of a vertical 
magnetic field, B z = B, and modes with purely vertical wavevectors, k z = k. We also take 
the rotation rate, pressure, and density to be constant on cylinders, i.e., Q(R), P(R), and 
p{R). We thus employ the same equations as those in BH02 (which are a simplification of 
the more general analysis in BH91). Although the problem we analyze is identical to that 
considered by BH02, we differ in the interpretation of the results. In particular, we identify 



-4- 



a set of modes that have all the features of growing convective modes in hydrodynamics, 
including inward transport of angular momentum. 

The thermal stratification of the flow can be described by the Brunt- Vaisala frequency, 

3 dPd\nPp-^ 3 n . 



N 2 = -N 2 R 



hp dR dR 



In the absence of rotation, the system is convectively unstable if N 2 > 0. In a rotating 
medium, but ignoring magnetic fields, convection is present only if iV 2 > k 2 , where k 2 = 
4f2 2 + dfl 2 /d In R is the epicyclic frequency.^ Note that N 2 > k 2 requires a sound speed 
comparable to the rotation speed. 

For a weakly magnetized medium, BH91 derived the behavior of linear perturbations 
including the effects of thermal stratification. For fluid displacements of the form £ oc 
exp[ikz + , ~ft], and for the particular geometry considered here, the growth rate 7 is given by 

1 



7 



-(kv A ) + 



N 2 -k 2 ± J(N 



K 2 ) 2 



16Q 2 {kv A ) 2 
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where v A = B / y/Anp is the Alfven speed. For the remainder of this section we focus on the 
growing mode in equation (§); this has 7 2 > and corresponds to the positive sign of the 
square root term. In the Appendix we present a more general analysis of both growing and 
oscillatory modes that helps clarify some of the issues discussed in this section. 

To understand the physics behind the dispersion relation in equation @ consider the 
equations for the radial and azimuthal displacements of fluid elements (BH91): 

/ dfl 2 N 



dt 2 



Solving equation (|) for 
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dt 2 ' — dt 

we can substitute into equation (§) to find 



^6 
dt 2 



where the three acceleration terms are of the form 

a H = (N 2 - k 2 )£r, a M , s = -(kv A ) 2 £ R , 



Q>ma 



m 2 k 2 v\ 
7 2 + k 2 v\ 



(3) 
(4) 

(5) 

(6) 



2 The requirements for convection in an unmagnetized rotating medium are the H0iland criteria (e.g., 
Tassoul 1978). These reduce to the condition N 2 > k 2 for the case considered here. For simplicity, we will 
refer to this condition as the H0iland criterion in this section. §2.1 treats the general case. 



- 5 - 



Since <9 2 £#/<9 2 t = 7 2 £_r it is straightforward to solve equation (|5|) and obtain the dispersion 
relation in equation (g). 

The first acceleration term in equation (|5|), a#, is a pure "hydrodynamic" or "H0iland" 
term (hence the subscript H) which is present even in the absence of magnetic fields. This 
term can be either stabilizing or destabilizing; it is stabilizing when N 2 < k 2 (the medium 
is convectively stable by the H0iland criteria) and destabilizing when N 2 > k 2 (convectively 
unstable). The second and third terms in equation (|6[) owe their existence to the magnetic 
field (hence subscript M). One of them is always stabilizing (subscript s) and the other is 
always destabilizing (subscript d). 

Equations ©-© have a very clean physical interpretation, as we explain in the rest of 
this section. For the simple geometry considered here (vertical field and vertical wavevector), 
gas and magnetic pressure forces are unimportant. If iV 2 = 0, the only forces in the problem 
are magnetic tension and rotating frame dynamics (e.g. the Coriolis force). Equations (^) 
and (|J) then describe the behavior of fluid elements that are coupled by the tension of 
magnetic field lines. As Balbus & Hawley (1992, 1998) have discussed, these equations are 
identical to those describing fluid elements coupled by a spring of spring constant (Icva) 2 - For 
a stratified medium with N 2 ^ 0, the radial acceleration of a fluid element has an additional 
component due to the radial buoyancy force. This is the N 2 £r term in equations (H) and (H). 
When N 2 < k 2 , this term does not have a large effect and introduces modest changes to the 
values of numerical factors. However, when N 2 > k 2 , which corresponds to a convectively 
unstable system by the H0iland criterion, the acceleration term an becomes positive and 
destabilizing. This introduces qualitatively new effects. 

The angular momentum flux of a linear mode is given by RQTr^, where Tr^ is the R — cj) 
component of the fluid stress tensor, 

Tr^ = p{8v R 8v 4> - 8v AR 8v A <p), (7) 

<5v is the perturbed fluid velocity, and <5v A = 5B/ ' y/Airp. Using the eigenfunctions in BH91 
we confirm BH02's expression for the stress tensor, namely 



p\5v 2 R \ D7 



k 2 v\ 
7 2 



dlnQ 



dhiR 



K 2 

20? 




where D = 1 + k 2 v\/'~f 2 and is a useful dimensionless stress whose sign gives the sign of 



We are now in a position to elucidate the physics of the linear instability calculation. 
Figure la shows a plot of the growth rate as a function of the dimensionless wavevector 
kvA/Q for 5 values of the Brunt- Vaisala frequency: N 2 /Q 2 = 0,0.5, 1, 1.5, and 2. We have 



- 6- 



assumed a Keplerian rotation profile, for which k 2 = Q 2 . Thus, the first two choices of iV 2 
correspond to a stable entropy gradient by the H0iland criterion, the third is neutrally stable, 
and the last two choices correspond to an unstable entropy gradient. Figure lb shows the 
dimensionless stress tj^ as a function of kvA/Q for the same five values of A 2 . It is useful 
to note that the self-similar ADAF model gives N 2 /Q 2 = 15 (7^ — l)/47 arf , where 7^ is the 
adiabatic index of the gas (Narayan & Yi 1994). | For example, if 7^ = 1.5, N 2 /Q 2 = 1.25. 

For A 2 = 0, the results reproduce the usual MRI. The growth rate is very small for 
kvA *C Q and peaks when kvA is ~ fi, with a peak value of 7 = 3f2/4; the angular momentum 
flux is always outwards (Tr^ > 0). The instability is clearly triggered by the magnetic field 
via the destabilizing term a,M,d, since the other two terms in equation (|^) are both stabilizing 
(negative). The MRI survives with qualitatively the same behavior even for non-zero A 2 , so 
long as the H0iland criterion for convective stability (iV 2 < k 2 ) is satisfied. 

Consider, however, a flow with iV 2 > k 2 . Such a flow is unstable to convection according 
to hydrodynamics, and is therefore the case of interest for CDAFs. Figure 1 shows that, 
for A 2 > k 2 , modes with kvA ~ are not very different from their A 2 = counterparts. 
The angular momentum flux is still outwards and the growth rate of the mode is comparable 
(though somewhat larger). By contrast, long- wavelength modes are very different. As Figure 
1 shows, long-wavelength modes are strongly unstable and they transport angular momentum 
inwards. In fact, as can be seen by inspection, the kvA *C limit of the dispersion relation 
in equation (with the + root) gives 7 2 = A 2 — k 2 , and the corresponding stress tensor 
is Tjtyj \p\5vr\ 2 ) = —k 2 /[2VI(N 2 — k 2 ) 1 / 2 ] < 0. The long-wavelength modes are clearly 
independent of the magnetic field and behave very differently from MRI modes. 

The physics of the kvA "C modes is as follows. In the pure MRI problem, namely 
A 2 = (or more generally A 2 < k 2 ), long- wavelength modes are only weakly unstable 
because the magnetic tension forces that are central to the MRI are weak; the field lines are 
hardly bent by a long- wavelength perturbation. By contrast, when A 2 > k 2 , the fluid has 
two destabilizing forces, the buoyancy term (a#) due to an unstable entropy gradient and 
the standard MRI term (a,M,d)- The buoyancy force that drives convection is independent of 
the wavelength of the perturbation, whereas the MRI depends on k. For long wavelengths, 
buoyancy is much more important and controls all the physics of the mode, both the growth 
rate and the angular momentum transport. Moreover, we show in the Appendix that the 
"convective" mode driven by buoyancy is the only unstable long-wavelength mode in MHD 
for a H0iland-unstable medium (A 2 > k 2 ). In particular, long- wavelength perturbations 
for which magnetic tension dominates (referred to as — modes in the Appendix) are stable 



3 This expression is not valid if j a d = 5/3 since then f2 2 = in the analytical models. 
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oscillatory waves rather than unstable modes (see eq. X3 and Fig. 3). 

Figure 2 explicitly compares the relative importance of the two destabilizing terms, the 
buoyancy term an and the magnetic tension term ciM,d- Let us identify a mode as an "MRP 
mode if an ^ a>M,d and as a "convective" mode if an > a>M,d- With this natural identification, 
we see that, for N 2 < k 2 , all unstable modes are MRI modes, regardless of the wavelength. 
Even when N 2 > k 2 , short- wavelength modes with kvA ^ ^ are still MRI modes. However, 
unstable long-wavelength modes with kvA <C ^ are clearly convective modes. The growth 
of these modes is due entirely to the unstable buoyancy force; correspondingly, the modes 
transport angular momentum inwards (Fig. lb), just as indicated by hydrodynamic studies 
(e.g., Ryu & Goodman 1992). Indeed these modes are virtually indistinguishable from their 
hydrodynamic counterparts. All of the properties of hydrodynamic convective instabilities 
therefore carry over to a magnetohydrodynamic analysis. 



2.1. The Hydrodynamic Limit 



Magnetic tension should be weak for fluctuations with kv a "C ^ regardless of the details 
of the magnetic field geometry and thermal stratification. In this subsection, we show this 
explicitly by considering the long-wavelength limit of the most general axisymmetric MHD 
dispersion relation. 

Balbus (1995) considered the linear, adiabatic, and axisymmetric stability of a dif- 
ferentially rotating weakly magnetized plasma in the presence of both vertical and radial 
stratification. He showed that the dispersion relation is given by 



/■••-• 



cJ 2 



^(DP)(D\n[Pp-^}) + ±- 3 D(RW) 



4fi 2 (k- v A ) 2 = 0, 



(9) 



where 



D 



(10) 



'knd_ _ _d_\ 
J^dz aR/' 

and Cj 2 = — 7 2 — (k ■ v A ) 2 (for fluctuations oc exp[7t + ik ■ r]). Unlike in the previous 
section, Q, P, and p are now allowed to be functions of both z and R, the magnetic field 
has an arbitrary direction (i.e., Br, B^, and B z components), and the wavevector is given 
by k = kfjTl + k z z. 



For k • v A = (but k/k z , k R /k z , etc., finite), equation 



1 



f-(DP)(D\n[Pp^}) + jr 3 D(R^ 2 



reduces to 
= 0. 
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Equation ([11]) is independent of the magnetic field. Indeed, as noted by Balbus (1995), 
it describes the linear, adiabatic, axisymmetric, hydrodynamic stability of a differentially 
rotating and thermally stratified medium (e.g., Goldreich & Shubert 1967). As is readily 
confirmed, equation ( ) has unstable solutions if and only if the H0iland criteria are satisfied 
(e.g., Tassoul 1978). 

The hydrodynamic limit in equation (|TT| ) formally requires setting k • va = in the 
MHD dispersion relation (eq. ||). If k • v A ^ 0, the full dispersion relation (eq. ||), 
which has twice the number of modes (because it is fourth order in 7 rather than second 
order), must be used. However, if the medium is unstable by the H0iland criteria, and if 
we focus on long-wavelength MHD instabilities with k • va <C fi, then the MHD system will 
behave like its hydrodynamic counterpart; in particular, the only unstable mode is effec- 
tively hydrodynamical and its growth rate is given by equation QTTJ) with small corrections 
~ (k • va/O) 2 <C 1. Note that equation fllTf ) describes the unstable growing mode only when 
the medium is H0iland unstable. If the medium is H0iland stable, the only possible insta- 
bilities are intrinsically MHD in nature and equation (^[) must be used even for k ■ va <C fl. 
This point is explained in more detail in the Appendix. 

The analysis in this section shows that all of the properties of axisymmetric hydrody- 
namic instabilities carry over to MHD in the long- wavelength limit. In particular, the results 
of the previous section are general, and are not an artifact of the simplifying assumptions 
made there. Although our quantitative analysis is restricted to axisymmetry, we suspect 
that long-wavelength non-axisymmetric modes will behave hydrodynamically as well (since 
they also do not significantly bend the magnetic field lines). 

3. Implications for CDAFs 

To assess the implications of our linear analysis it is useful to consider wavelengths 
relative to the scale height of the disk, H pa c s /Q, where c s is the sound speed. The 
maximally growing MRI mode has kH ~ P 1 ^ 2 , where j3 ~ c 2 s /v\ is the ratio of the thermal 
energy to the magnetic energy in the disk. By contrast, the longer wavelength convective 
modes have kH <C /3 1//2 ; more quantitatively, the linear analysis suggests that buoyancy 
forces become important for kH ^ 0.3/3 1 / 2 (Fig. 2). For these modes to be of interest they 
must fit in the disk and so must have kH ^ 1. This in turn requires > 1, i.e., that the 
magnetic field must be sub-equipartition (/3 ^ 10 may suffice). 

Given the above analysis, we propose the following identification between the hydro- 
dynamic CDAF model and the present MHD results. If the magnetic field saturates at a 
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sub-equipartition value, the most unstable MRI mode operates on a scale ~ f3~ l l 2 H <C H 
and transports angular momentum outwards. Moreover, if the medium is H0iland unsta- 
ble, only such short-wavelength fluctuations transport angular momentum outwards. Longer 
wavelength convective fluctuations are dominated by scales ~ H and transport angular mo- 
mentum inwards (Fig. 1). 

Which fluctuations are more important? Or, more precisely, which part of the power 
spectrum of fluctuations is more important, those with kH ~ 1 or those with kH ~ f3~ 1 ^ 2 l 
This cannot really be addressed analytically because we do not understand the saturation of 
the linear instabilities (and it is clear that they are coupled; e.g., convection, by stirring the 
fluid, can help to build up the magnetic field). Nonetheless, it is typically the case that the 
energy in a turbulent plasma is dominated by the largest unstable scales in the medium. This 
suggests that, if (3 3> 1 (perhaps (3 k, 10 is sufficient), the larger-scale convective fluctuations 
will dominate the dynamics, as proposed in the hydrodynamic models. Moreover, the inward 
angular momentum transport by convective modes depends on the magnitude of N 2 — k 2 
(Fig. lb). Thus there is a natural way for the inward convective transport to self-adjust 
to counteract outward transport by small scale MRI fluctuations. In the limit of efficient 
convection and (3 3> 1 we would expect marginal stability to convection according to the 
H0iland criteria, just as in the hydrodynamic models (Quataert & Gruzinov 2000). 

It is important to stress several points: 

(1) Although we believe that our linear stability analysis sheds important light on the 
physics of radiatively inefficient accretion flows, it does not directly address the fully de- 
veloped nonlinear turbulence that is of primary importance to the accretion flow structure. 
This requires numerical simulations. 

(2) The picture proposed here is very similar to the hydrodynamic models of CDAFs; the 
primary change is that we have theoretical support from a magnetohydrodynamic analysis. 

(3) Balbus (2001) has shown that magnetized dilute plasmas can be convectively unsta- 
ble in the presence of an outwardly decreasing temperature, rather than just an outwardly 
decreasing entropy. It would be interesting to incorporate this into future work. 

(4) The extent to which one can usefully distinguish scales in the power spectrum that 
are convective from those that are magnetorotational depends on (3, and becomes more and 
more useful if (3 3> 1. It would be straightforward if the convective and MRI modes occupied 
distinct and well-separated regions of wavevector space. Then even the nonlinear physics of 
the instabilities could be distinct. To cite an example from plasma physics, heat transport 
in fusion devices is due to instabilities on both proton and electron Larmor radii scales and 
one can often distinguish the physics of each mechanism separately even in the nonlinear 
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regime (e.g., Dorland et al. 2000; Rogers, Dorland & Kotschenreuther 2000). As Figure 1 
shows, such a separation in wavevector space is not present for our problem. On the other 
hand, Figs. 1 and 2 do show quite clearly that the physics of the convective modes is very 
different from that of the MRI modes. 

(5) For a H0iland-unstable medium (iV 2 > k 2 ), there is a nice continuity between 
the hydrodynamic limit and the strong field MHD limit. In wavevector space, the field- 
affected modes (the MRI modes) are restricted to larger values of k where magnetic tension 
is important. For very weak fields, these modes are at very large k, and most of the long- 
wavelength modes that probably dominate the dynamics are convective. As the field strength 
increases, however, the MRI modes expand downwards in k space and become more and 
more important. Another way to look at this is as follows. When the fluid is unstable 
according to the H0iland criterion, the switch that determines whether or not convection is 
important is not the presence or absence of a magnetic field, but rather the magnitude of the 
dimensionless parameter kvA/Q ~ kH/fi 1 / 2 . When kvA/Q is small, convection dominates 
the physics, whereas when it is ^ 1, the magnetic field dominates. 



4. Conclusion 

In this paper we have used linear stability theory to assess the conditions under which 
a hydrodynamic analysis is appropriate for an intrinsically magnetohydrodynamic accre- 
tion flow. We have shown that, when the medium is unstable by the H0iland criteria, 
long-wavelength instabilities with kvA <C ^ are effectively hydrodynamical. Physically, this 
is because magnetic tension is unimportant for sufficiently long- wavelength perturbations. 
Formally, the dispersion relation for a magnetohydrodynamic plasma reduces in this limit to 
that of hydrodynamics (§2). This shows that all of the properties of hydrodynamic instabil- 
ities carry over to an MHD analysis in the long- wavelength limit. 

For the particular case of radiatively inefficient accretion flows that are unstable to 
both convection and the MRI, the relevant parameter that distinguishes convective behavior 
from magnetorotational behavior is the quantity fciu/^ ~ H / '(A/3 1 / 2 ). Long- wavelength 
modes for which X/H >> (3~ 1 ^ 2 are convective, while modes for which X/H ~ j3~ x l 2 are 
magnetorotational. The convective modes are primarily driven by buoyancy, not magnetic 
tension (Fig. 2). Moreover, they transport angular momentum inwards rather than outwards 
(Fig. lb), just as indicated by hydrodynamic studies. For the long-wavelength modes to 
be of interest the accretion flow must have (3 ^> 1 so that the convective fluctuations fit 
inside the disk. Under these conditions we believe that the hydrodynamic CDAF models 
can accurately describe the structure of the accretion flow. The precise (5 that allows a 
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hydrodynamical analysis is uncertain and can be determined only by magnetohydrodynamic 
simulations (the linear analysis suggests that (3 ^ 10 may suffice). 

It is interesting to note that the saturation (3 of the MRI does not appear to be uni- 
versal. Given the plausible dependence on (3 highlighted above, this may have important 
implications for the structure of radiatively inefficient accretion flows. Three-dimensional 
MHD simulations find that that if the field is initially toroidal a radiatively inefficient ac- 
cretion flow develops into a CDAF-like state (Igumenshchev, Abramowicz & Narayan 2002), 
while if the field has a vertical component the structure is very different (Hawley, Balbus & 
Stone 2001; Hawley & Balbus 2002; see also Stone & Pringle's 2001 2D simulations). 

To conclude, it is important to emphasize that the primary conclusion of the analytic 
arguments in this paper is not that CDAFs actually do describe the structure of radiatively 
inefficient accretion flows, but rather that they could. In particular, our linear MHD analysis 
supports all of the key CDAF assumptions so long as (3 ^> 1. 

RN, IVI, and EQ were supported in part by NSF grant AST-9820686, RFBR grant 
00-02-16135, and NASA grant NAG 5-12043, respectively. 
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A. Appendix 



In the main text we focused exclusively on unstable modes and therefore considered 
only the positive sign of the root in equation (0). Here we generalize the discussion and 
consider both unstable and stable modes. To this end, we rewrite equation (|2|) as follows 



7 2 = -(kv A ) 2 + ± '- 



1 ± Jl + 16n 2 (kv A ) 2 /(N 2 - k 



2\2 



(Al) 



The + and — signs in equation flA~I] ) map directly onto the + and — signs in equation 
(0) when A^ 2 — k 2 > 0, but the mapping is reversed when A^ 2 — k 2 < 0. Equation ([TO]) 



is, of course, completely equivalent to equation (g), but for this Appendix we find it more 
convenient to work with the ± convention in equation ( |A1| ) because of its behavior when 
A^ 2 — k 2 changes sign. In what follows, we refer to the + sign in equation ([Al]) as the "+ 
mode" and the — sign as the "— mode." 

Let us focus on the long-wavelength limit kv A /Q <C 1. The growth rates of the + and 
— modes are then given by 

7 2 - (iV 2 - n 2 ) + - (kv A ) 2 , (A2) 




7* « - + 1 ) (kv A f. (A3) 

In the limit kv A /Q <C 1, the growth rate of the + mode is the same as in hydrodynamics: 
7+ = A^ 2 — k 2 . If A^ 2 > k 2 (H0iland unstable), 7 2 is positive and the mode is unstable, while 
if A^ 2 < k 2 , 7 2 is negative and the mode is oscillatory in nature. In either case, the mode is 
completely insensitive to the magnetic field. 

The — mode, on the other hand, has a growth rate that vanishes as kv A — > 0, which is 
a clear indication that this mode is strongly influenced by the magnetic field. Indeed, the 
mode owes its very existence to the presence of the magnetic field. The — mode is associated 
with the classic MRI of BH91. 

An inspection of equations (|A2|) and (|A3|) reveals the following interesting result: as- 
suming that |A^ 2 — k 2 \ < 4f2 2 (which is almost always the case), only one of the two modes 
can be unstable. When the medium is convectively stable (A^ 2 < k 2 ), the + mode (which 
we argued above is hydrodynamic in nature) is obviously stable. In this case, the — mode 
is unstable and gives the long-wavelength limit of the MRI. However, when the medium is 
convectively unstable (A^ 2 > k 2 ), then it is the + mode that grows, while the — mode is 



-In- 



stable. Thus, in the long wavelength limit, the MRI actually disappears if the medium is 
unstable by the H0iland criterion. Fortunately, this is precisely when the medium is con- 
vectively unstable so there are always unstable modes available, regardless of the sign of 
N 2 - k 2 . 

Figure (|3|) shows 7 2 as a function of {kvA) 2 for the + and — modes for representative 
values of N 2 /Q 2 . We see that modes that are stable (7 2 < 0) in the limit kvA — ► remain 
stable for all values of k, while modes that are unstable at small values of kvA reach a 
maximum growth rate for kvA ~ Q and then become stable with increasing kvA- 
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Fig. 1. — (a) Dimensionless growth rate as a function of dimensionless wavevector kvA/Q 
in a Keplerian system (k, 2 = Q 2 ), for five choices of the Brunt- Vaisala frequency: N 2 /Q 2 = 
(solid line), 0.5 (dotted line), 1 (short-dashed line), 1.5 (long-dashed line), 2 (dot-dashed 
line). The system is unstable by the H0iland criterion for N 2 /VL 2 = 1.5 and 2. For these 
two cases, the growth rate remains finite in the long- wavelength limit kvA/Q <C 1 and the 
modes are indistinguishable from the corresponding unstable hydrodynamical modes, (b) 
Dimensionless shear stress corresponding to the same modes as in (a). For N 2 /Q 2 = 1.5, 
2, and fcf^/fi 1, the stress is negative, implying that these convective modes transfer 
angular momentum inwards, just as in hydrodynamics. 
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Fig. 2. — Variation of a H /aM,d versus dimensionless wavevector fcw^/fi for N 2 /Q 2 = 
(solid line), 0.5 (dotted line), 1.5 (long-dashed line), 2 (dot-dashed line). The ratio a,H/Q>M,d 
measures the importance of the hydrodynamic buoyancy force relative to the leading desta- 
bilizing force due to magnetic fields (see eqs. For the two cases in which the system 
is convectively unstable according to the H0iland criterion, viz., N 2 /Q 2 = 1.5, 2, the buoy- 
ancy force clearly dominates at long wavelengths, confirming that these modes are primarily 
convective, not magnetorotational. 
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Fig. 3. — The dimensionless growth rate 7 2 /^ 2 versus dimensionless wavevector kv^/Q for 
N 2 /Q 2 = 0, 0.5, 1.5, 2. A mode with a positive value of 7 2 is unstable, while one with a 
negative value of 7 2 is stable (oscillatory). The solid and dashed lines correspond to the + 
and — modes, respectively, in equation fl^l]). Each curve is labeled by the value of N 2 /fl 2 . 



